Radial flow has little effect on clusterization at intermediate 
energies in the framework of the Lattice Gas Model 
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The Lattice Gas Model was extended to incorporate the effect of radial flow. Con- 
. trary to popular belief, radial flow has little effect on the clusterization process in 

intermediate energy heavy-ion collisions except adding an ordered motion to the 
particles in the fragmentation source. We compared the results from the lattice gas 
£^ ■ model with and without radial flow to experimental data. We found that charge 

, yields from central collisions are not significantly affected by inclusion of any rea- 

t> 

' sonable radial flow. 
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PACS numbers: 25.70.Pq, 24.10.Pa, 64.60.My 

Assuming local equilibrium and ignoring any collective motions, many statistical models 
are used to explain the cluster formation process in multifragmentations. Ignoring certain 
collective motions, such as sideward flow, is often justified by the short range nature of 
correlations which dominates the cluster formation process. Sideward flow only adds a 
uniform boost to the momentum of the particles in the neighborhood of a given particle and 
thus does not change the phase space distributions in the local rest frame. However, essential 
difference arises when the collective motion adds a nonuniform boost to the particles in the 
the neighborhood of a given particle. Such nonuniform boost could not be cancelled by a 
simple change of reference frame, and could induce a nontrivial change in the phase space 

n 

distribution. One of the simple nonuniform boosts is radial flow an d we will consider 
the effects of it on the cluster production process. 

Experimentally, we know there is a strong radial flow in central heavy-ion collision. The 
radial flow could often account for half of the kinetic energy of the particles. In the paper , 
essentially different pattern of charge yields are found in central collisions than in peripheral 
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collisions. Kunde et al. [2j compared two systems with similar thermal energy (at least the 
authors of claim so), but quite different radial flow energy. One of the systems is from 
low energy central collisions, where significant radial flow is present. Another system is from 
higher energy peripheral collisions where little radial flow is present. The charge yields from 
the central collisions are found to be of exponential type while those from the peripheral 
collisions are of power law type. Whether such a strong effect for radial flow exists or not 
has become a key issue for the understanding of cluster formations in intermediate energy 
heavy-ion collisions. 
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and the thermodynamic canonical model 
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12L ll3L Il4j , the radial flow 



can be added only a posteriori. Thus, the radial flow is not included at a fundamental level. 



Transport models, such as BUU and QMD could incorporate t 
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re flow features naturally, 

n 

19j . In this sense, 
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but have trouble to reproduce the details of charge yields 
the lattice gas model bears some advantages in being able to incorporate flow in a natural 
way and also simulates some fine details of cluster productions. 

In a previous paper, we have developed the lattice gas model to include radial flow j^. 
It was also applied to a small system with mass A = 84 to test the effect of radial flow. 
The radial flow incorporated in the calculation is quite small, on the order of 2 AMeV, and 
the effect of radial flow was found to be very small. In this paper, we will test the radial 
flow effects in much larger systems where the flow effect is presumably much larger, and 
also the experimental data are available Q]. And we will show that a moderate radial flow 
will have little effect on the cluster formation process in the current model. To attribute the 
different patterns found in central collisions to radial flow motion, we need to have radial 
flow energies much larger than the total energy of the collision system. 

For completeness of discussion, we will briefly review the microcanonical lattice gas model 
[^l, and the incorporation of radial flow in the model (see also Then we will try to 

reproduce the experimental data from central and peripheral collisions in the lattice gas 
model without any flow. We also present the calculations with quite large radial flow. A 
summary of the essential finding in this work is given in the end. 

The idea of the lattice gas model is to discretize space and put the nucleons on the 



lattice 



22L |23| . Each lattice has at most one nucleon to simulate the short range repulsion 



effect between nucleons and also partially Pauli principle. The attractive NN interaction is 
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simulated with a bond between a pair of neighboring nucleons, where the bond strengths 
are take as e np = —5.33 MeV, e nn = e pp = 0. These bond values could roughly reproduce 
the average binding energy of clusters and the energy density of infinite nuclear matter. The 
phase space factor is difficult to calculate directly, and is usually sampled with a Metropolis 
Algorithm. We are using the microcanonical ensemble so the total energy E to t is constant. 
We start with a given configuration, hence a given E pot . Thus the kinetic energy for a given 
E tot is Ek in = E tot — E pot . The available phase space of A = N + Z nucleons have this kinetic 
energy is known analytically: 

^ -Efy n = fiiM (2m)3A/2 ^ • (1) 

We call this value the constant energy integral Q(Ekin)- We now attempt a switch in 
configuration space. The potential energy will change to E' t . To conserve energy, the 
kinetic energy of the system should be E' kin = E tot — E' pot . If Q(E kin )/Q(E kin ) > 1, the 
switch is accepted. If it is less than 1, the probability of accepting such a switch is given by 
the ratio. Since all configurations have equal weights, this satisfies the principle of detailed 
balance. After many such switching attempts, an event is accepted. We now have to assign 
momenta to the nucleons so that the total kinetic energy of the A nucleons add up to correct 
rinetic energy. This can be done by standard Monte-Carlo sampling (for further details, see 



The inclusion of radial flow is quite straight forward in lattice gas model. Radial flow 
is a collective motion and does not change the thermal motion. So we may add the flow 
motion at the end of configuration sampling. Each nucleon acquires additional momentum 
Pf'(i) = c[f(i) — fcm), and the constant c is adjusted to give the required flow energy. With 
radial flow, the total energy is now E tot = E therm + E pot + Ef [ow . From this, it seems that 
flow has only made the part of the total excitation energy unavailable for thermal motion. 
But this is not totally true. The flow motion also affect the cluster formation in the current 
model. 

After the sampling of configurations and momentum assignment, we can search for clus- 
ters. Two neighboring nucleons are considered to belong to the same cluster when the energy 
for this pair is lower than zero: 

E pair = £ + e<0, (2) 
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where p r is the relative momentum and e is the bond between the pair. This clusterization 



prescription is standard in Lattice Gas Model |2 



the Coniglio-Klein prescription 
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and was shown to be similar to 



241 ] . A cluster formed under this condition is, by 



definition, stable against particle emissions and can be compared with the charge yields from 



heavy- ion reactions |22j, |2£| . The particle stable property of this clusterization method was 
also verified in a molecular dynamical simulation [2y]. 

In the pair definition, Eq.(J2J), the relative motion has components from radial flow, so 
that in general, pairs are less bound than without radial flow. This effect tends to reduce the 
cluster size formed on the lattice. On the other hand, increasing thermal motion also has 
the effect of reducing the binding between neighboring nucleons and subsequently reduces 
the cluster size. So the question we have to ask is, qualitatively how big is the radial flow 
effect relative to the thermal effect. 
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FIG. 1: The charge yields from experiment and from lattice gas model for two different systems. 
The data from experiment (j^]) are represented by full symbols, while the results from lattice gas 
model are shown with open symbols. The two systems used in our simulations are, one with A = 
150, Z = 60, N = 90, E therm = 7.5 AMeV on 8 3 lattice and one with A = 320, Z = 130, N = 190, 
Etherm = 13 AMeV on 10 3 lattice respectively. The charge yields for the A = 150 system show 
a typical power law behavior and that for the A = 320 system show a exponential behavior. All 
charge yields are normalized to the charge yield of Li element. 
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Here we will present the results from microcanonical lattice gas model in two large sys- 
tems, for which the experimental data were reported before |2j. The two systems we used 
are: A = 150, Z = 60, N = 90 on 8 3 lattice which corresponds to peripheral Au + Au 
collisions at Ei a b = 1 GeV/nucleon; A = 320, Z = 130, iV = 190 on 10 3 lattice which 
corresponds to central Au + Au collisions at Eiab = 100 MeV/nucleon. First, we will try 
to reproduce the experimental data in the lattice gas model with only thermal motion and 
no radial flow. As shown in FigQ the full symbols are data for the two system while the 
open symbols are the corresponding simulation results. Overall, the lattice gas model can 
reasonably produce the general features of the data. 

For Lattice gas results on A = 150 system with E t herm = 7.5 AMeV thermal energy, 
the charge yields show a power law behavior, and closely follow the trend of data except 
yields for the highest charges. We find the lattice gas model overestimates the charge yields 
when the cluster size is approaching half of the system size, Z = 30; this maybe due to edge 
effect. For the A = 320 system, the lattice gas results show a nice exponential behavior for 
charge yields at a higher excitation energy E t herm = 13 AMeV. Even though the energy scale 
correspondence in the lattice gas model does not correspond exactly to that in a real nuclear 
system, we find these thermal energies are still in reasonable agreement with expectations. 
For peripheral reaction as used in the experiment (A = 150 system), typical excitation 
energy is estimated to be around 6 — 8 MeV (see |2| and reference therein). For central 
collisions, typical excitation energy is expected to be around 15 MeV. 

Following the argument of |2j, we have tried to find a reasonable fit to the data on 
A = 320 system by fixing the thermal energy to 7.5 AMeV while changing the radial flow 
energy. But no matter what reasonable flow energies we used the charge yields are almost 
not changed within statistical error in the simulations. This is a surprising result. As we will 
see, random thermal motions are order of magnitude more efficient in breaking the bonds 
between neighboring nucleons than radial flow motion. 

The effect on charge yields from radial flow and thermal motion are shown in FigEl For 
the charge yields from the lattice gas model, the difference between 12 and 13 AMeV pure 
thermal motion is very large. If we want to produce the charge yields of E t h erm = 13 AMeV 
by introduce a large radial flow on the E therm = 12 AMeV results, we need about 72 AMeV 
radial flow. In other words, the random thermal motion is 72 times more efficient than radial 
flow motion at breaking the bonds between neighboring nucleons. 
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FIG. 2: The charge yield for the system A = 320, Z = 130, N = 190 on 10 3 lattice in the 
microcanonical lattice gas model. Filled symbols are for simulations with only thermal motion and 
no radial flow at all, at Etherm = 12 and 13 AMeV respectively. Open symbols are for simulations 
with thermal energy fixed at Etherm = 12 AMeV while radial flow energy is varied between 24, 48, 
and 72 AMeV. As can be seen, the charge yields from simulation with Etherm + Efi ow = 12 + 72 
AMeV roughly match the charge yields with pure thermal motion E t herm = 13 AMeV. All charge 
yields are normalized to that of Li element. 

The reason flow affects clusterization so little is because bond is between nearest neigh- 
bors. For a bond to break we must have large (pi — pjj). With flow we add an extra 
momentum difference of c(f% — r 2 ) = cl where |/| = (6.25Y^fm. This is rather small and 
does not depend upon how big the dissociating system is. We will now show that Etherm 
is much more efficient in breaking bonds than Efi ow . We do this by first calculating the 
average value of (pi — p 2 ) 2 with pure E therm . We then consider if instead of E therm , the same 
energy was put into Efi ow . We find that the value of (pi — p 2 ) 2 falls by about a factor of 50. 

For thermal motion ((pi — P2) 2 ) = 2(pf) = 2{p 2 ). It is reasonable that 

(f) = 2m^p . (3) 
where A is the number of nucleons in the system. The validity of this equation can be 
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established for our model by turning the discrete problem into a continuous one. Then 

Etherm = « / 4:71— p 2 dp, 

A = k 4np 2 dp, 
Jo 

J P p 2 dp 



2, .2 



We now want to estimate momentum mismatch when E t h eTm = E^ ow = £\ Turn 
this discrete problem again into a continuous problem. We assume the dissociating system 
is uniformly distributed in 3 times the normal volume (A more realistic density distribution 
is given in j27|). Calling this expanded radius R and the reduced density p we find 

c 2 A R 5 

Etherm — Eflow ~ ^ 4lX p — 

2m 5 

This leads to c 2 = l0mE t herm/3AR 2 . Comparing 2 x 2mE t herm/A with c 2 l 2 , we find the 
former is larger by a factor 45 if the dissociating system has 320 nucleons. A more elaborate 
estimate using the geometry of the lattice will give a factor of 50. 

From our simulation, we see this factor is 72 instead of 50, an increase of 44%. In fact, 
the above estimate only gives a lower limit for this relative efficiency, any deviation from a 
uniform distribution of nucleons in the simulation will definitely increase this ratio. Further, 
the additional thermal energy for the pair of neighboring nucleons is actually weighted with 
a distribution for the pair of neighboring particles, so that resulting average will be higher 
than our estimation. The overall increase of 44% in the ratio of relative bond-breaking 
efficiency is most likely the result of a combination of these factors. 

Since the flow values we used here are a lot higher than the total energy deposition in 
central collisions at E\ ah = 100 MeV/nucleon, the lattice gas simulations suggests that there 
is little effect of radial flow on the clusterization process at intermediate energies. The main 
effect that caused the charge yields to behave like exponential type rather than power law 
type is the increase of thermal energy in central collisions. 

The inefficiency of radial flow on breaking clusters is also observed by Gulminelli and 
Chomaz 23| i n a canonical lattice gas model. For other commonly used thermal models, 
which have the same assumption of sudden dissociation, a simple inclusion of radial flow as 
we did here will not likely produce a significant difference. Different assumptions about the 
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cluster formation may change the relative efficiency of the radial flow and thermal motion 
in breaking up clusters, but the large difference between radial flow and thermal motion will 
hold in any reasonable models. 

The above discussion did not address the issue of the dynamical evolution of the heavy- 
ion reaction system, especially during the violent stage of the reaction. The clusterization 
method used in Lattice Gas Model is similar to the r — p phase space coalescence. If we 
blindly apply the estimate to normal density, we would find that thermal motion is 22 times 
more efficient in breaking bonds than radial flow motion. Thus we may expect our idea to 
shed some light on the violent stage of the reaction as well. The dynamical models has also 
made significant progress to in reproducing cluster formations. In QMD model, algorithms 
like ECRA and SACA have been formulated to identify asymptotic clusters at an early 
stage |3, [lfl 0|. In BUU model, stochastic Langevin method was introduced to address 
the stochastic nature of cluster formations These methods are calculational intensive. 
We hope future dynamical simulations will confirm our expectations. 

To summarize, we have tested the effect of radial flow on cluster formation in microcanon- 
ical lattice gas model. We used quite large system where the radial flow effect is expected to 
be strong and data already existed. We found that reasonable amount of radial flow could 
not significantly change charge yields in the simulation. For the two systems as discussed in 
Q], we found the charge yields changed from power law behavior to an exponential behavior 
purely because of increased thermal energy. Thermal energy is order of magnitude more 
efficient in reducing the size of large clusters than radial flow. These conclusions should be 
valid for most thermal models where dynamical evolution of the source is ignored. 
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